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Abstract 


I review standard big bang nucleosynthesis and some versions of non¬ 
standard BBN. The abundances of the primordial isotopes D, He-3, and 
Li-7 produced in standard BBN can be calculated as a function of the 
baryon density with an accuracy of about 10%. For He-4 the accuracy is 
better than 1%. The calculated abundances agree fairly well with obser¬ 
vations, but the baryon density of the universe cannot be determined with 
high precision. Possibilities for nonstandard BBN include inhomogeneous 
and antimatter BBN and nonzero neutrino chemical potentials. 

1 Introduction 

Big bang nucleosynthesis (BBN) is among the main observational evidence for 
big bang. The discovery of the cosmic microwave background (CMB) provided 
us with the temperature scale of the early universe, and allowed the calculation 
of the primordial nuclear abundances produced in the big bang. The four light 
isotopes, D, 3 He, 4 He, and 'Li are mainly produced in the big bang, and the 
calculated abundances agree fairly well with astronomical observations. 

Standard big bang nucleosynthesis (SBBN) has a single free parameter, the 
baryon-to-photon ratio, 



(1) 


which is related to the present baryonic contribution to the critical density fib 
via the Hubble constant Hq = /il00 km/s/Mpc by 


77io = 10 1O 77 = 274 Q b h 2 . 


( 2 ) 


For decades, BBN has provided the best determination of the amount of bary¬ 
onic matter in the universe. The agreement with observations is obtained in the 
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range 7710 = 1.5 ... 6 . Despite optimistic claims from time to time, BBN has not 
really progressed towards a much more precise determination of 77 . Observers 
claim higher precision from determinations of primordial abundances of single 
isotopes, but disagree with each other or, within the context of SBBN, with 
primordial abundances of other isotopes. Difficult questions about systematic 
errors in observations and chemical evolution relating the present abundances 
to primordial abundances have prevented further progress. 

During the past year, a competing method for estimating the amount of 
baryonic matter has appeared. In the angular power spectrum of the anisotropy 
of CMB, the relative heights of the even and odd acoustic peaks are sensitive to 
the baryon-to-photon ratio. The first preliminary estimates from the Boomerang 
g§ and Maxima-1 Q experiments appeared to be in conflict with BBN, giving 
a higher baryon density, flbh 2 ~ 0.03, or 7710 ~ 8 j|(J. The Boomerang collab¬ 
oration has since revised their estimate downward, to f hh 2 = 0 . 022 ±q qo 3 J38|, 
which agrees with SBBN, but the Maxima-1 estimate has been revised upward 
to Qbh 2 = 0.0325±0.0125 (95% c.l.) J53) . With the coming satellite experiments 
CMB may surpass BBN as the method for estimating 77 . BBN will then become 
a tool for understanding the astrophysics of chemical evolution, by telling us 
the primordial abundances. 

While SBBN is simple and natural, and is at present in reasonable agreement 
with observations, there is interest in studying nonstandard variants of BBN. 
For one thing, BBN is a sensitive probe of the physics of the early universe. 
If we change something about our assumptions regarding the conditions in the 
early universe, or the physics relevant at that time, we are likely to change 
the primordial abundances and ruin the agreement with observations. Thus for 
many things BBN provides the strongest constraint. 

On the other hand, from time to time there have been suggestions for dis¬ 
agreement between the estimated primordial abundances of the different iso¬ 
topes, and/or other ways of estimating 77 . If such disagreements persist, non¬ 
standard BBN (NSBBN) may be the solution. 

I shall begin with a review of the physics of SBBN, and then discuss a few 
NSBBN scenarios. 

2 Physics of Big Bang Nucleosynthesis 

In the early universe the temperature is falling as the universe expands. The 
time scale depends on the number of particle species which are relativistic at 
that time. In the standard case these are electrons, positrons, photons, and 3 
species of neutrinos. 

There is lots of radiation and very little matter. The amount of baryonic 
matter is not known exactly, and is given by 77 , the only free parameter in SBBN. 

Weak reactions are converting neutrons into protons. At first the reaction 
rate is high enough to maintain the equilibrium neutron-to-proton ratio 



( 3 ) 
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which is falling with temperature. Here £ e = l^u e /T is the electron neutrino 
degeneracy parameter (p Ve is the electron neutrino chemical potential). For 
£ <C 1, we have roughly 


£ ~ 


n v - n p 
n 7 


(4) 


In SBBN we assume homogeneous conditions with £ e ~ 0. 

At a temperature of about T ~ 0.8 MeV the neutrinos decouple and after 
that the neutron abundance evolves via free neutron decay 


n —> p + e + v e . 


(5) 


Nuclear reactions begin by neutrons and protons producing deuterium. This 
reaction keeps the deuterium abundance close to its equilibrium value. Because 
of the large amount of photons in the background radiation and the low bind¬ 
ing energy, Bd = 2.22 MeV, of deuterium, the deuterium abundance does not 
become large until the temperature has fallen to about 70 keV. Only then can 
the nuclear reactions proceed beyond deuterium. 

As the temperature falls further the Coulomb barrier shuts down the nuclear 
reactions. Because of the short time and bottlenecks due to the lack of stable 
nuclei at masses A = 5 and A = 8, the reactions do not proceed beyond A = 7. 

There are 10 important reactions which take nucleosynthesis beyond deu¬ 
terium: 


d(p, 7 ) 3 He d(d,p)t d(d,n ) 3 He 3 H e(n,p)t 

t(d, n) 4 He 3 He(d,p) 4 He 3 He(a, 7 ) 7 Bc t(a, 7 ) 7 Li 

7 Be(?r,p) 7 Li 7 Li(p, a) 4 He 

Afterwards tritium /3-decays into 3 He and 'Be becomes 7 Li by electron capture. 

Since 4 He has the highest binding energy per nucleon (for A < 7), almost 
all neutrons end up incorporated in 4 He. Thus the primordial 4 He abundance 
Y p is determined by the neutron fraction at nucleosynthesis time. This in turn 
is determined by the competition between the weak reaction rates and the rate 
at which the temperature falls due to the expansion of the universe. A higher 
baryon density causes nucleosynthesis to take place at a higher temperature, 
when there are more neutrons left. Thus Y p increases with ry. 

A higher ry means more efficient nuclear burning to 4 He, leaving less of 
the “impurities”, D and 3 He. There is also less directly produced 7 Li, but a 
higher density allows more 7 Bc to be produced, so that the final primordial 7 Li 
abundance has a more complicated dependence on ry (see figure 1). 

3 Accuracy of Nucleosynthesis Yields 

The accuracy of the SBBN calculation of produced primordial abundances de¬ 
pends on how accurately the various thermonuclear reaction rates are known. 
These rates cannot be calculated from first principles at present, and are thus 
based on laboratory measurements. BBN codes make use of compilations of 
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these rates @ 0001 - The Caughlan&Fowler || 1988 compilation was 
updated on some key BBN rates by Smith et al. Q. The latest compilation is 
by the NACRE collaboration Q. 

For BBN, the most significant difference between the new NACRE rates and 
the older rates is that deuterium production is now higher. With the NACRE 
rates, the O’Meara et al. ]39j result D/H = 3.0 ± 0.4 x 10 ~ 5 corresponds to 
ryio = 5.9 ± 0.5 instead of ryi 0 = 5.6 ± 0.5 obtained with the older rates. 

Based on how accurately the reaction rates are known, Cyburt et al. |J have 
estimated the accuracy of SBBN yields in the range 7710 = 1... 10 to be better 
than 0.3 % for Y p , 3-7 % for D, 3-10 % for 3 He, and 12-19 % for 'Li. 


4 Observations 


I shall leave the more detailed discussion of the observations relating to pri¬ 
mordial abundances to other participants of this workshop, and just list various 
recent results in the literature. 

For 4 He, Olive et al. [Q combined results of different observers and esti¬ 
mated a primordial abundance Y p = 0.234 ± 0.003. Izotov&Thuan ]h| used 
their own large set of observations to arrive at Y p = 0.244 ± 0.002. The errors 
are 1 -a statistical errors. There has been much recent discussion of possible 
systematic errors. Peimbert et al. jidj have recently measured the 4 He abun¬ 
dance in the SMC, whose proximity allows a better control of systematic errors. 
They obtained Y = 0.2405 ± 0.0018 for the SMC, and from this they estimate 
Y p = 0.2345 ± 0.0026 for the primordial abundance. Thuan&Izotov |55| have 
now refined their estimate to Y p = 0.2443 ± 0.0015. After discussing various 
systematic effects, they conclude that because of systematic errors this could be 
an underestimate by ~ 2-4%. 

The observed deuterium abundance in the local interstellar medium is D/H = 
1.50 ±0.10 x 10~ 5 Q, but there is evidence for local variations J 57 J. For exam¬ 
ple, Sonneborn et al. Q have recently reported D/H = 2.18±g;fj r x 10 -5 along 
one line of sight. Since deuterium is destroyed in stellar processes, its primordial 
abundance should be greater than the local abundance. 

The most promising method of obtaining the primordial deuterium abun¬ 
dance is the measurement of D/H from Lyman-series absorption by high-redshift 
clouds. Based on three such measurements and upper limits from others, 
O’Meara et al.jdfj obtain D/H = 3.0 ± 0.4 x 10 -5 . Measurements on one 
such system suggest a higher deuterium abundance, possibly larger than 10~ 4 
[58, 56 j. Recently, the deuterium in two more such systems has been observed 
[ 10, 32, ^2), supporting the low D/H of O’Meara et al. p9| . 

The estimates for the primordial 7 Li abundance are based on the rather 
uniform abundances in population II halo stars. Bonifacio&Molaro ]|| ob¬ 
tained log 10 ( 7 Li/H) = —9.80 ± 0.012 ± 0.05 and Ryan et al. |47j obtained 
log 10 ( 7 Li/H) = —9.88 for the mean abundance. This measured abundance is 
close to the minimum 7 Li from SBBN (see figure 1). It may be depleted some¬ 
what from the primordial abundance by stellar processing. Pinsonneault et 
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Figure 1: The SBBN yields of the light isotopes as a function of the baryon- 
to-photon ratio rj. The rectangles corresponds to the various observational esti¬ 
mates for the primordial abundances mentioned in the text: 4 He is from Olive 
et al. (4^] and Izotov&Thuan JltJ, D is from O’Meara et al. ]3!| and 7 Li is from 
Ryan et al. Q. For 7 Li a conservative upper limit log 10 ( 7 Li/H) < —9.4 is also 
shown. Figure from [^0[. 
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al. Q estimate log 10 ( 7 Li/H) = —9.65 ... — 9.25 for the primordial abundance. 
The observed ' Li may also include a galactic contribution, so that the primor¬ 
dial abundance could be even lower. For the primordial abundance Ryan et 
al. [Q estimate log 10 ( 7 Li/H) = —10.04... —9.72 and Suzuki et al. Q estimate 
logicn Li/H) = -9.97 ■ • • - 9.77. 

Comparing these estimates for primordial abundances to the SBBN yields 
(figure 1) we see that there is some tension between D/H, which favors a higher 
baryon density, and Y p and 7 Li/H, which favor a lower baryon density. The best 
estimate from CMB now agrees with the high ?? from D/H. Thuan&Izotov ^5| 
now conclude that their value for Y p is in good agreement with this when one 
allows for the systematic error. If one accepts the Thuan&Izotov J55| result for 
Y p the remaining disagreement is with the 7 Li. Possibly 'Li processing is not 
yet understood well enough, and the tighter limits above are too stringent. 

5 Nonstandard BBN 

There are many proposed possible modifications to SBBN. I shall go over just 
a few of these NSBBN scenarios. 

A higher energy density at a given temperature would lead to faster expan¬ 
sion in the early universe. This would be caused by the presence of additional 
relativistic particle species, e.g., additional neutrino species. The shorter time 
scale would mean that there are more neutrons left at nucleosynthesis time, and 
thus more 4 He is produced. This is the most important effect. Since a higher 
Y p tends to lead to worse agreement with observations, BBN sets an upper limit 
to the speed-up possible. This upper limit is usually parameterized in terms 
of the effective number of (light) neutrino species N v . Different authors get 
different upper limits depending on the observational constraints chosen. Two 
recent results are N„ < 3.2 by Buries et al. Q and N v < 4 by Lisi et al. Q. 

The effect on other isotopes is roughly that their abundance curves vs. the 
baryon density are shifted towards higher r]. The higher density compensates 
for the shorter time available. This is a smaller effect than the effect on 4 He. 

A large neutrino degeneracy, £ ^ 0, would increase the energy density in 
neutrinos at a given temperature, also leading to a speed-up of the expansion 
rate. If the degeneracy is in electron neutrinos, a much larger effect is that 
( e /0 shifts the balance of weak reactions (see eq. ||). If £ e > 0, so that there 
are more electron neutrinos than electron antineutrinos, we get less neutrons 
and thus a lower Y p . Since this has a large effect on Y p but a small effect on the 
other isotopes, we can use £ e to “dial in” a desired value of Y p , as £ e is otherwise 
unobservable. 

One can combine the above scenarios by having both a significantly faster 
expansion rate and £ e > 0, so that the effects on Y p cancel each other to keep it 
in the observational range, while the effect of the speed-up on the other isotopes 
remains pf| . If one wants to use neutrino degeneracy only to achieve this speed¬ 
up, this scenario requires |£ M | £ e or |£ T | £ e . This way larger ?? are allowed 

than in SBBN. However, since one is relying on a small effect, the large effects 
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Figure 2: Results from IBBN for f l^h 2 = 0.030 (7710 = 8.22). This is the density 
suggested by the preliminary Boomerang and Maxima-1 results, r is the distance 
scale giving the separation between centers of high- and low-density regions, and 
f v is the volume fraction of the high-density regions. The distance scale is given 
in comoving units at T = 1 MeV. (1 km at 1 MeV corresponds to 1.9 x ICC 4 pc 
at present.) The thick lines show the observational constraints Y p < 0.248 (to 
the left of the curve) and D/H = 2.9-4.0 x 10 -5 . IBBN can thus bring this 
high 77 into agreement with 4 He and D observations. Thin lines are contours 
of log 10 ( 7 Li/H), showing that 7 Li remains problematically high. This figure is 
for a “spherical shell” geometry which is more successful in reducing the 7 Li 
yield than other geometries tried. The dotted curve shows where figure 3 cuts 
through this figure. Figure from J2f|. 


having cancelled each other, this scenario requires a very large speed-up for 
a significant effect. A large N„ leads to other cosmological effects which put 
limits to this scenario. In particular, the recent CMB anisotropy results place 
an upper bound to N v [|l5| , |l 6 | . When the CMB constraints are combined with 
this BBN scenario [^2[ gdf gif |l^] the conclusion is that the upper limit to 7710 
can be raised from the SBBN 7710 < 6 to 7710 < 7... 8 only (24J . 


6 Inhomogeneous BBN 

Another way to modify SBBN is to relax the homogeneity assumption, and 
allow 77 to be a function of position at nucleosynthesis time. In the usual kind of 
inhomogeneous BBN (IBBN), the baryon density is positive, 77 > 0, everywhere, 
but one can also consider the possibility of having antimatter regions, where 
77 < 0. This latter case I call antimatter BBN (ABBN), and I discuss it in the 
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Figure 3: Like figure 2, but now we show the dependence on the third IBBN 
parameter, the density contrast R between the high- and low-density regions. 
The dotted line shows where figure 2 cuts through this one. The solid thick line 
is the constraint Y p < 0.248 (above the line) and the dashed thick lines give the 
D/H constraint. The thin lines are the log 10 ( 7 Li/H) contours. This model can 
just reach below our “conservative upper limit” log 10 ( 7 Li/H) < —9.4. Figure 
from 












next section. 

The crucial question in IBBN is the distance scale r of the inhomogeneity. 
If the mechanism causing the inhomogeneity is connected with inflation, inho¬ 
mogeneity at any scale could be produced. The observed isotropy of the CMB, 
however, requires 77 to be highly homogeneous at distance scales larger than 
about 10 Mpc. Thus IBBN cannot be used to explain, e.g., different deuterium 
abundances at different high-redshift Lyman-absorbers. 

In the usual IBBN scenarios one assumes a much smaller distance scale of 
the inhomogeneity, so that while 77 is inhomogeneous at nucleosynthesis time, 
the matter from different regions gets mixed later, resulting in a homogeneous 
baryon density with homogeneous abundances, which however are different from 
the SBBN abundances for the same 77 . 

In the simplest version of IBBN, each small region undergoes SBBN with 
its own baryon-to-photon ratio, and one gets the final abundances by averaging 
over the SBBN abundances over the distribution in 77 . Leonard&Scherrer [fiof 
have shown that this way one can relax, or even remove, the observational lower 
limit to 77 , but the upper limit to 77 remains unchanged. 

If the mechanism producing the inhomogeneity was not connected with in¬ 
flation, causality requires the distance scale to be smaller than the horizon at 
that time. A favorite candidate for producing the inhomogeneity has been the 
QCD phase transition at T ~ 150 MeV, when the horizon was about 1 pc (co¬ 
moving). The order of the phase transition is not known, but if it was of first 
order, it would have proceeded through nucleation of bubbles of hadronic matter 
in the ambient quark-gluon plasma. Near the end of the transition there would 
have been shrinking droplets of quark-gluon plasma, where much of the baryon 
number would have been concentrated due to the difficulty of baryon number 
crossing the phase boundary. This way very high density contrasts could be 
produced. The relevant distance scale is close to the neutron diffusion scale at 
nucleosynthesis time. An inhomogeneity of this scale would produce a strongly 
inhomogeneous neutron-to-proton ratio due to neutron diffusion out of the high- 
density regions, leading to a drastic change in BBN n HI- In the first papers 
the effect was overestimated; for more accurate calculations see, e.g., 0 P~ill - 

A more extensive list of references is given in the review articles by Reeves [|I4| 
and Malaney&Mathews Q and in, e.g., Q. For a range of distance scales one 
can get less 4 He and more D than in SBBN, favoring a higher 77 . One can also 
get some reduction in the ' Li yield, but not much. Because of the large number 
of parameters in the scenario, quantitative results are not easily summarized. 
The IBBN scenario is contrasted with recent observations in |2^] and |2f| . 

The properties of the QCD phase transition can be estimated by lattice QCD 
calculations. While a distance scale interesting for IBBN can not be ruled out, 
unfortunately a too small distance scale appears more likely. 
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Figure 4: Antimatter regions. Figure from 



7 BBN with Antimatter Regions 

We do not know the origin of the baryoasymmetry, rj > 0, of the universe. 
Theories of baryogenesis try to explain this excess of matter over antimatter and 
to predict (postdict) the observed value of rj. While the simplest models usually 
lead to a homogeneous rj, in many baryogenesis models the baryon density 
may come out inhomogeneous, and in some regions the asymmetry may have 
the opposite sign, so that after local annihilation we end up with regions of 
antimatter 0 - 

Annihilation will then proceed at the matter-antimatter boundary. The 
smaller the antimatter regions the sooner they are completely annihilated. From 
the spectrum of the CMB we know that there was no major annihilation going 
on close to recombination time, and the observed cosmic diffuse gamma ray 
radiation puts tight limits on annihilation after recombination. This leaves 
us three possibilities not in contradiction with observations: 1) the antimatter 
regions annihilated well before recombination time, 2) the amount of antimatter 
was much less (a factor of about 1(W 6 or less) than the amount of matter, or 
3) the antimatter regions are separated by large distances, comparable to the 
present horizon, or larger. 

The first possibility leads to an interesting NSBBN scenario, antimatter 
BBN (ABBN) |h| |4h], jl6|, [2(], ^?], ^9) . We consider antimatter regions which 

are larger than 10^“ pc but smaller than 100 pc (comoving). These annihilate 
after T = 1 MeV, but before recombination, and the annihilation process can 
significantly affect nucleosynthesis, or modify the abundances afterwards. Since 
in the end we have complete annihilation of the antimatter, we need to assume 
an excess of matter over antimatter, so that the antimatter/matter ratio R < 1. 

Antimatter regions smaller than 10 -2 5 pc annihilate before nucleosynthe¬ 
sis, when the temperature of the universe is T = 1 MeV- 70 keV. The mixing 
of matter and antimatter is due to (anti)neutron diffusion, and therefore the 
annihilation is preferably on neutrons. This reduces the 4 He production. The 
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Figure 5: Observational constraints for ABBN with 7710 = 6 . R is the antimat¬ 
ter/matter ratio, and ta is the radius of the antimatter region, given in comoving 
units at T = 1 keV. (1 nr at 1 keV is 4.26 x 10 6 m = 1.38 x 10 ~ 10 pc today.) 
The shaded region is the region allowed by the constraints Y p = 0.228-0.248 
and D/H = 2.2-4.0 x 10 -5 . ABBN is thus able to remove the tension between 
D and 4 He observations. This region is also allowed by other constraints: the 
spectral distortion of the CMB and upper limits to 3 He/D; the excluded regions 
are above and to the right of these curves. Because 6 Li is very fragile, it is 
difficult to make definitive conclusions about its primordial abundance from ob¬ 
servations. Therefore we just show two contours, 6 Li/ 7 Li = 0.1 and 6 Li/ 7 Li = 1. 
The 7 Li yield is fairly constant, close to the SBBN value 7 Li/H = 2.8 x 10 -10 , 
over the whole figure. Figure from p9[. 


observational lower limit to Y p therefore provides an upper limit to the amount 
of antimatter R. 

Larger antimatter regions survive until nucleosynthesis, which consumes the 
free neutrons. This stops further annihilation for a while. There is then a 
second burst of annihilation well after nucleosynthesis, when proton diffusion 
becomes effective, at T < 3 keV. The annihilation process then changes the 
nuclear abundances through several mechanisms. 

Gamma rays from annihilation photodisintegrate nuclei. The major effect is 
the production of D and especially 3 He from photodisintegration of 4 He. 

Antiproton annihilation on helium nuclei produces 3 H, 3 He, D, protons and 
neutrons. These nuclear remnants have high energies, and the energetic 3 H lead 
to nonthermal production of 6 Li by the endoergic reaction t(a , ?r) 6 Li. 

Thus D, 3 He, and 6 Li abundances are increased over the SBBN yields. There 
is no big effect on 7 Li. 
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Figure 6: Same as the previous figure, but with rjio = 8. ABBN can accom¬ 
modate higher values of r) than SBBN, except for the 7 Li constraint, which is 
essentially the same for ABBN and SBBN. Figure from [ [f9| . 


At T < 1 keV, the energy released in annihilation does not get fully thermal- 
ized, but results in a distortion of the CMB spectrum. Since no such distortion 
has been observed, CMB places constraints an the amount of antimatter an¬ 
nihilating below T ~ 1 keV, but before recombination. Near T ~ 1 keV the 
universe is still strongly radiation dominated, so the energy release in annihi¬ 
lation is small compared to the energy in the background radiation. As the 
temperature falls, the matter-to-radiation energy ratio increases. Therefore the 
limits from CMB to the amount of antimatter become progressively stronger for 
larger distance scales, for which the annihilation occurs later. For scales smaller 
than 1 pc, we get, however, a stronger limit from 3 He overproduction. 

Since ABBN reduces Y p and increases D/H, it allows for a larger baryon 
density, at least from those observational constraints. The constraint from 7 Li, 
however, remains essentially unchanged. 


8 Conclusions 

Standard BBN is compelling in its simplicity. While there is controversy among 
the observers and some apparent discrepancy between the estimated primordial 
abundances of the different isotopes and SBBN, these are probably not serious, 
and most likely represent difficulties in making observations and estimating 
primordial abundances based on observed ones. SBBN is thus not in trouble. 
Unfortunately, because of these difficulties, it is not able to pin down the baryon- 
to-photon ratio very precisely. It is somewhere in the range 7710 = 1.5-6, or 


12 





flbh 2 = 0.005-0.022. The high redshift deuterium measurements point towards 
the upper end of this range, 7710 ~ 5-6, or Ll^h 2 ~ 0.02. 

The recent estimates from CMB anisotropy, f l^h 2 = 0.022+qqq 3 from Boom¬ 
erang |3f| and il^h 2 = 0.0325 ± 0.0125 (95% c.l.) from Maxima-1 53 also favor 
this upper end of the range. We are eagerly waiting for more precise CMB 
measurements in the coming years. 

While standard BBN is in good shape, there is interest in studying non¬ 
standard BBN: to assess the robustness of SBBN, to constrain possibilities for 
nonstandard physics or cosmology, and to be ready to provide relief if observa¬ 
tional discrepancies turn out to be serious for SBBN. 

We discussed here four NSBBN scenarios: 1) electron neutrino degeneracy, 
2) electron neutrino degeneracy combined with a speed-up of the expansion 
rate due to extra energy density, 3) inhomogeneous BBN, and 4) antimatter 
BBN. All these scenarios are able to relieve the tension between the D and 
4 He observations. The three last ones may also allow a larger baryon density 
than SBBN, but with some difficulty: 2) is constrained in that respect by other 
cosmological constraints and 3) and 4) cannot do much for the 7 Li constraint. 

I thank Elina Sihvola for permission to reproduce figures from pH, ptj. 
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